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The results of a lattice simulation of time-separated charge overlap for the charged pion are discussed. The 
expected result exp[—{Eq — mTy)t] for large charge density overlap time separations, t, is clearly visible in 
the Fourier transform, indicating that the elastic limit can be achieved at low to medium momentum values on 
present-sized lattices. The implications of this result for direct lattice simulations of hadron structure functions 
are discussed and a brief presentation of the lattice formalism is given. 



1. INTRODUCTION 

In trying to better understand the internal 
structure of hadrons in lattice QCD, charge over- 
lap techniques are playing an increasingly im- 
portant role. These methods are based upon 
the simulation of the hadronic matrix elements 
< /i(O)|r[j;J(r,i)J^(0)]|/i(O) >, where J;J^" are 
d,u flavor current densities. In their original 
form[|l],|j these measurements, using t = be- 
tween the currents, allowed qualitative studies of 
lattice hadrons to be carried out. However, it 
has been pointed out that large Euclidean time 
separations of the currents allows one to extract 
form factor data from this matrix element in the 
so-called elastic limit Q. This is interesting in it's 
own right, but it is also clear that by replacing the 
flavor currents in this matrix element with the full 
electromagnetic ones, we are then studying the 
basic matrix element needed for hadron structure 
functions0|. Thus, in moving from the qualita- 
tive to the quantitative stage, the elastic limit of 
these matrix elements forms an important bridge, 
leading to both elastic and inelastic properties of 
hadrons. The question of how readily this limit 
may be implemented in lattice simulations is the 
subject of the present investigation. 



where 

^ < ^+(O)|T[V(r + x,tK(x,0)]k+(O) > . 

X 

(2) 

pU,d g^j,g |-jjg ^ flavor charge densities and 'T 'is 
time-ordering. Figure 1 is a symbolic represen- 
tation of the measurement in Eq. (^). Assuming 
t > and inserting a complete set of states, this 
results in 



X 



<X(q)K(0)|7r+(O)>e 



(3) 



For large Euclidean times the sum reduces to a 
single term and we find in the continuum limit 
that 



Qoo(q^^) 



t»i 



(4) 



when the SU(2) flavor symmetry is present. 



2. THEORY 

Consider the lattice Euclidean time-separated 
charge overlap distribution for zero momentum 
pions: 

Qot(q',i)^Ee"^'''"^oo(r,t), (1) 



Figure 1. Symbolic representation of the d, it 
charge overlap measurement. 



Thus, by separating the currents in time it is in 
principle possible to damp out the intermediate 
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state contributions and to measure form factors 
at arbitrary lattice momenta. On a finite space- 
time lattice, however, achieving this elastic limit 
is problematical. The exponential damping fac- 
tor, Q-(Eq-m^)t ^ jg j^Q^ large and the fixed time 
positions of the initial and final pion interpolating 
fields limit the possible charge density separation 
times. We now turn to a numerical investigation 
of this limit. 

3. SIMULATION 

The present discussion is based on the numer- 
ical results found in Ref. where more calcu- 
lational details may be found. This study was 
conducted on 12 quenched /3 — 6.0 SU(3) con- 
figurations {16^ X 24) at K = .154 (the largest 
K value studied in Ref. ||].) The pion interpo- 
lating fields were located at time slices 4 and 21 
and the current densities were positioned as sym- 
metrically possible in time between these sources. 
Two Wilson quark propagators per configura- 
tion, with origins at the positions of the inter- 
polating sources, were necessary to extract the 
relative overlap function V^Qir^t) in Eq. (||). 
(7^QQ'(r, i) is actually the large Eucldean time 
limit of a similarly denoted quantity in Ref. [Q . ) 
The pion interpolating field at time slice 21 was 
constructed with quark propagators which were 
smeared over the entire spatial volume using the 
lattice Coulomb gauge. Both point and smeared 
pion fields were used at time slice 4. These propa- 
gators produce four-point functions which project 
exactly on zero pion momentum. Charge density 
self-contractions were neglected in forming these 
quantities. 

Fig. 2 represents a log^o plot of the 
Fourier transform of these distribution func- 
tions, QQo{<l^,t), at the two lowest lattice spa- 
tial momentum values, |q| = tt/8 (upper points), 
^/2{n/8) (lower points), as a function of rela- 
tive time separation between current densities. 
The solid lines shown in this figure come from 
the expected asymptotic exponential falloff spec- 
ified by Eq. (||), using the (smeared) k = .154 
data from Table 1 (m^r — .369, mp = .46) of 
Ref. Q , assuming the vector dominance form for 
the pion form factor: F^(g^) = 1/(1 + g^/m^). 



Figure 2. Logj^Q plot of the Fourier transform, 
Qqq, of the spatial density charge correlation 
function as a function of relative time separation, 
t, between charge density operators. 



(We also assume the continuum relation Eg = 
(toJ -I- q^)^/^.) Actually shown in this fig- 
ure are results for both point-to-smeared (O) as 
well as smeared-to-smeared (□) correlation func- 
tions. In all cases, the expected functional depen- 
dence ~ g-(s<j-m„)t jg pj-gggut, indicating that 
by time step 7 or 8 single exponential behavior 
has emerged. This is similar to the number of 
time steps needed in hadron spectrum calcula- 
tions. This behavior is remarkable because al- 
though we are damping out intermediate states 
as t increases, we are also moving closer to pos- 
sible contaminations from the fixed interpolating 
fields at either time end. In fact, we do not see 
any indications of such contamination in the data. 
These correlation functions are also unusual be- 
cause the asymptotic approach is from below, in- 
dicating damping of negative terms in the d, u 
correlation function. Ref. ||] indicates that in the 
continuum limit, these are primarily positive G- 
parity states (t). 

Although the point-to-smeared and smeared- 
to-smeared results exhibit essentially overlapping 
error bars, the smeared-to-smeared values are 
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Figure 3. Extracted pion form factor, F^^, as a 
function of q'^ /rrr^. The solid line is vector domi- 
nance. 



systematically low compared to the point-to- 
smeared results, indicating a slight dependence 
on the form of the interpolating field used. How- 
ever, the results of Ref. j|] indicate that such a 
dependence decreases that is, as the 

physical regime is approached. 

By fitting time steps 7 through 10 of Fig. 2 with 
a single exponential of known slope and removing 
the kinematical factor in Eq. (Q), a value of FT^{q'^) 
can be obtained. The results for the two calcu- 
lated /rn^p values using point-to-smeared corre- 
lation functions are shown in Fig. 3 (O). The val- 
ues found are consistent within errors with vector 
dominance, shown as the solid line. Also shown 
in this figure are the results (□) from a previous 
thrce-point-function simulation of the pion form 
factor Comparison shows that the error bars 
of these different techniques are of the same or- 
der of magnitude for similar numbers of configu- 
rations. (Note the different (3 values and lattice 
sizes of the two simulations, however.) 



4. STRUCTURE FUNCTIONS 

In order to begin to understand the significance 
of these results for structure functions, let us re- 
call the basic form of these quantities (continuum 
Minkowski expression): 

W^p{q\v)^^ I d^xe-^'i-^ 
2tt J 

•^(7r+(O)|[J„(x),J^(0)]k+(O)). (5) 

We adopt the lab frame where the external pion 
momentums are zero and use the full electromag- 
netic currents. In this frame v = —p ■ q/m — go- 
In order to make contact with lattice expressions, 
we imagine discretizing space, but keeping the 
time variable continuous. Using the correspon- 
dences, 

|7r+(p)) ^ [Nsa^2Ej,]i\7r+{p)>, 

Tcont. —3 jlatt. 

where N s is the number of space points in the 
lattice, this gives (all quantities on the lattice) 

N S'°° 

• <7r+(O)|[J„(x),J^(0)]|7r+(O) > . (6) 

By a standard set of manipulations, this can be 
shown to result in the expression: 

W^p{q\,y) = N^J2^i:^ -Ex+ m^) 

X 

■ < 7r+(0)| J40)|X(q) X X(q)| J^(0)|7r+(O) > . 

(7) 

How is this quantity to be measured on the 
lattice? Let us consider the generalization of 
Eqs. (|l|) and (^ to arbitrary components of the 
full electromagnetic current. Define in Euclidean 
space 

r,x 

• <7r+(O)|r[J„(x + r,i)j0(x,O)]|7r+(O) > . (8) 
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Again assuming t > 0, this results in 

X 

■ < 7r+(0)| J„(0)|X(q) X X(q)| J;3(0)|7r+(O) > . 

(9) 

We notice that the connection between Eq. (^) 
and (0) is formally given by the inverse Laplace 
transform, 

-J pC~\~'lOQ 

Wo.p{q^v)^— dte'^'Qo.^iq^t), (10) 

where c > and the contour is closed in the 
left-half complex t-plane. Since we are at fixed 
spatial momentum q, the lattice data extracted 
will actually cut a parabola-like path given by 
X = /{2my^ q^ — q^) in the x, plane (a; = 
q^/{2mv) and = q^ — q^). In the lab frame 
the range of v when < a; < 1 is |q| < v < 
{Eq — m^) . One can assure these kinematical con- 
straints in the lattice model by restricting the as- 
sumed forms of QQ/3(q^, t) to assure contributions 
to Waf} (q^ , v) only within the continuum range 
of V. We also need to make a quasi-continuum 
assumption about the lattice data: quantities 
like Woo {'f , ^) should be very poorly fit by a 
positively weighted sum of exponentials. This 
assumption implies that Qoo (q^ , t) is given by 
products of exponentials and inverse powers of 
t. (In this case Eq. (^ becomes purely formal 
but Eq. (0) continues to hold.) This whole dis- 
cussion has been carried out for the pion, but it 
is clear that there is no barrier to applying these 
techniques to the phenomenologically more inter- 
esting proton case as well. 

5. DISCUSSION 

There are now two workable techniques for ex- 
tracting form factor data from lattice simulations: 
direct current insertion j|] and elastic charge over- 
lap. The elastic limit also makes it possible to 
perform direct simulations of hadron structure 
functions. Attention so far has been focused on 
the moments of such functions, which are given 
by the operator product expansion. These ex- 
pansions are based upon separation of the short- 
distance physics, calculated perturbatively, from 



the long-distance part, which can be evaluated 
on the lattice in the form of certain operator ex- 
pectation values. Because of the perturbative as- 
sumption these methods work best at large q^\ 
direct lattice simulations must use low dimension- 
less (qa)'^, so these two techniques should be com- 
plementary. The question of whether the direct 
method can reach the scaling regime is still open, 
but note that the chiral limit (elastic) q^ range in 
Ref. H was .6 < < l.QGeV^. Assuming the 
elastic limit of all necessary flavor-diagonal and 
non-diagonal four-point functions can be demon- 
strated, direct simulations of structure functions 
should be feasible with current computer technol- 
ogy. 
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